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O ABSTRACT 

(N 

Context. Because of their short evolutionary time-scales, the earliest stages of high-mass star formation prior to the existence of any 
O i embedded heating source have barely been characterized until today. 

^A Aims. We study the fragmentation and dynamical properties of a massive starless gas clump at the onset of high-mass star formation. 

Methods. Based on Herschel continuum data we identify a massive gas clump that remains far-infrared dark up to 100 [im wave- 
l/~) lengths. The fragmentation and dynamical properties are investigated by means of Plateau de Bure Interferometer and Nobeyama 

£SJ 45 m single-dish spectral line and continuum observations. 

Results. The massive gas reservoir (between ~800 and ~1600M o , depending on the assumed dust properties) fragments at spatial 

scales of ~ 1 8000 AU in four cores. Comparing the spatial extent of this high-mass region with intermediate- to low-mass starless cores 
•^T from the literature, we find that linear sizes do not vary significantly over the whole mass regime. However, the high-mass regions 

O squeeze much more gas into these similar volumes and hence have orders of magnitude larger densities. The fragmentation proper- 

ties of the presented low-to high-mass regions are consistent with gravitational instable Jeans fragmentation. Furthermore, we find 
multiple velocity components associated with the resolved cores. Recent radiative transfer hydrodynamic simulations of the dynamic 
collapse of massive gas clumps also result in multiple velocity components along the line of sight because of the clumpy structure of 
I the regions. This result is supported by a ratio between viral and total gas mass for the whole region <1. 

Conclusions. This apparently still starless high-mass gas clump exhibits clear signatures of early fragmentation and dynamic col- 
lapse prior to the formation of an embedded heating source. A comparison with regions of lower mass reveals that the linear size of 
star-forming regions does not necessarily have to vary much for different masses, however, the mass reservoirs and gas densities are 
orders of magnitude enhanced for high-mass regions compared to their lower-mass siblings. 

Key words. Stars: formation - Stars: early-type - Stars: individual: IRDC 18310-4 - Stars: massive - ISM: clouds - ISM: kinematics 
and dynamics 
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qq 1. Introduction strong fragmentation during earliest evolutionary stages (e.g., 

\Q Bontemps et al. 2010)? What are the kinematic properties of the 

. Independent of the various formation scenarios for high-mass gas? Are the clumps sub . or super _ v i ria i ? Do we see streaming 

stars that are discussed extensively in the literature (e.g., motions indicative of turbu lent flows (e.g., Bergin et al. 2004; 

O Zinnecker & Yorke 2007; Beuther et al. 2007), the initial con- Vazquez -Semadeni et al. 2006; Heitsch et al. 2008; Banerjee 

' ditions required to allow high-mass star formation at all are still t i 2009) 17 

poorly characterized. The initial debate even ranged around the _ , .... . , ... 

L .• u tu u- u .1 i . , . • . , To address these questions high spatial resolution is manda- 

>• question whether high-mass starless gas clumps should exist at „ , , , , 

. I_i ,, , .i ... „ e ■ , .. _. ■ tory. Furthermore, to study the dense gas, dust continuum as well 

all, or whether the collapse of massive gas clumps starts imme- J , ,. ' / & . ', _ 

r^ j- , i »u » 1*1 u ■ tu u- u ■ as spectral line observations are required. Therefore we tackle 

rN diately without any clear starless phase in the high-mass regime , F . . ... . ? r , • . • . 

'_h , \, .. . , onmx r, . ,. .■ • j- 7 tx. * a. Z- these questions via a dedicated study of the genuine high-mass 
~3 (e.g., Motte et al. 2007). Recent studies indicate that the time . M . mT . n ..,.„ . ' b _ . , * , , 

^3 , 6 • i ■ t tU t u aa a , f starless gas clump IRDC 18310-4 that was first identified by 

where massive gas clumps exist without embedded star forma- _, . „ & , ,„„„,.. _, . , , , , _ J 

,..,,.,., , (tmm n , . ., Sndharan et al. (2005). Their early study based on 1.2 mm con- 

tion is relatively short (on the order of 50000 yrs), but neverthe- . . . v . ' ..,,„, . . , , 

, , • , . , , , . ., „ i.i tinuum emission, 8 um extinction and NH3 emission revealed an 

less, high-mass starless gas clumps do exist (e.g., Russeil et al. . . ,,, ; •, ^ r.^«-.« 

oaia rf 1 u » 1 omiN infrared dark gas clump with a mass of ~ 840 M Q , a tempera- 
2010; Tackenberg et al. 2012). . . ° , ^ T¥T ,, ,. ,. ... °, . . F . 

_ . . d . ,, »■•, , ture of ~ 18K, and a NH3(1,1) line-width measured with the 

Current questions in that field are: Are high-mass gas clumps rlt , , lnn tl ..,.„„ c , _, -1 rp, , ■ 

, . , / . , . , • , , Ettelsberg 100 m telescope with 40 or ~ 1.7kms . The kine- 

dominated by a single fragment or do we witness already . ,. ° c ., . „ „ _ , . /T , , 

J ° ° J matic distance of the source is ~ 4.9 + 0.26 kpc (Ragan et al. 

— 2012a, recalculated following Reid et al. 2009) and its velocity 

Based on observations carried out with the IRAM Plateau de Bure e . oc«i_ -1 /c -ji. » 1 nnns\ 

mAw . ,, „.„„„.,„„ ,„ . .„„ of rest vi sr ~ 86.5kms (Sndharan et al., 2005). 

Interferometer. IRAM is supported by INSU/CNRS (France), MPG v ' 

(Germany) and IGN (Spain). The data are available in electronic form The target region is at a projected separation of ~ 3.4' (at 

at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or the given distance corresponding to ~ 5.1 pc) from the lumi- 

viahttp://cdsweb.u-st4rasbg.fr/cgi-bin/qcat?J/A+A/ nous high-mass protostellar object (HMPO) IRAS 18310-0825 



H. Beuther et al. : Fragmentation and dynamical collapse of the starless high-mass star-forming region IRDC 183 10-4. 

Fig. 1. Compilation of the continuum data from 70/im to 500 /im wavelength in color-scale as labeled in each panel. The scale is chosen in each 
image individually to highlight the most important features. Contour levels of the 870 yum data start at the 3<x levels of 0.189mJybeairT 1 and 
continue in in 3<x steps. The white small box in the top-left panel marks the smaller region discussed in the remaining figures. A scale-bar is shown 
in the bottom-left panel. 



(Sridharan et al., 2002; Beuther et al., 2002) but otherwise in 
no outstanding part of our Galaxy. This region is part of the 
Herschel EPoS key project (Earliest Phases of Star Formation) 
that studies a large sample of high- and low-mass star-forming 
regions at the earliest evolutionary stages (Ragan et al. 2012a, 
Launhardt et al. 2013, Henning et al. 2010, Linz et al. 2010, 
Beuther et al. 2010, 2012b, Stutz et al. 2010). 

Except for these characterizations, little additional informa- 
tion exist for that region. Here, we present Herschel far-infrared 
continuum observations as well as single-dish and interferomet- 
rically obtained N2H + (l-0) observations from the Nobeyama 
45 m telescope and the Plateau de Bure Interferometer. The latter 
data also allow us to derive high-spatial-resolution 3 mm contin- 
uum data. 



2. Observations 

2.1. Herschel 

The cloud complex with a size of ~ 6' x 6' was observed 
with PACS (Poglitsch et al., 2010) on Herschel (Pilbratt et al., 
2010) on April 19th, 2011 (PACS OBS ID 134221906[0-3]). 
Scan maps in two orthogonal directions with scan leg lengths 
of ~ 10' were obtained with the medium scan speed of 20"/s. 
The raw data have been reduced up to level- 1 with the HIPE 
software (Poglitsch et al., 2010; Ott, 2010), version 6.0, build 
1932. Besides the standard steps, we applied a 2 m 'level deglitch- 
ing, in order to remove bad data values from a given pixel map 
by cr-clipping the flux values which contribute to each pixel. We 
used the time-ordered option and applied a 25cr threshold. The 
final level-2 maps were processed using Scanamorphos version 
8 (Roussel, 2012). Since the field of view contains bright emis- 
sion on scales larger than the map, we applied the "galactic" 
option and included the non-zero-acceleration telescope turn- 
around data. The flux correction factors provided by the PACS 
ICC team were applied. The 70 //m extinction feature and the 
PdBI continuum and line data exhibited a small offset of ap- 
proximately 3" in R.A. and 1" in Declination. Although the in- 
terferometer positional accuracy is usually superior because of 
the extremely well-known positions of the reference quasars, 
we double-checked the 70 //m PACS data via a comparison to 
the Spitzer MIPSGAL 24 fim survey which has also a very 
high spatial accuracy (Carey et al., 2009). Again, the same off- 
set between MIPSGAL 24/j.m and Herschel 70/im was found. 
These Herschel data were taken at still an early Herschel time, 
and positional offsets on this order are known during that time. 
Therefore, we shifted the 70 //m map by this little offset in the 
high-resolution overlays from Figure 2 onwards. The spatial res- 
olution of the 70, 100 and 160 //m data is 5.6", 6.8" and 11.4", 
respectively. 

Maps at 250, 350, and 500 /um were obtained with SPIRE 
(Griffin et al., 2010) on March 11th, 2010 (SPIRE OBS ID 
1342192067). Two times two 20' scan legs were used to cover 
the source. The data were processed up to level- 1 within 
HIPE, developer build 5.0, branch 1892, calibration tree 5.1 
using the standard photometer script (POF5 .pipeline .py, dated 
2.3.2010) provided by the SPIRE ICC team. The resulting level- 
1 maps have been further reduced using Scanamorphos, version 



9 (patched, dated 08.03.2011). This version included again the 
essential de-striping for maps with less than 3 scan legs per scan. 
In addition, we used the "galactic" option and included the non- 
zero-acceleration telescope turn-around data. Corrections for the 
transition from point sources to extended sources were applied 
according to the SPIRE Manual. The spatial resolution of the 
250, 350 and 500 yum data is 18.1", 24.9" and 36.6", respec- 
tively. 



2.2. Plateau de Bure Interferometer 

We observed IRDC 18310-4 with the Plateau de Bure 
Interferometer during five nights in October and November 
2009 at 93 GHz in the C and D configurations covering pro- 
jected baselines between approximately 13 and 175 m. The ob- 
servations were conducted in a track-sharing mode together 
with the IRDC 18454-1 in the neighborhood of W43 (Beuther 
et al., 2012b). The 3 mm receivers were tuned to 92.835 GHz 
in the lower sideband covering the N2H + (l-0) as well as the 
3.23 mm continuum emission. At the given frequency, the pri- 
mary beam of the PdBI is ~ 54". For continuum measurements 
we placed six 320 MHz correlator units in the band, the spec- 
tral lines were excluded in averaging the units to produce the 
final continuum image. Temporal fluctuations of amplitude and 
phase were calibrated with frequent observations of the quasars 
1827+062 and 1829-106. The amplitude scale was derived from 
measurements of MWC349 and 3c454.3. We estimate the final 
flux accuracy to be correct to within ~ 15%. The phase refer- 
ence center is R.A. (J2000.0) 18:33:39.468 and Dec. (J2000.0) 
-08:21:07.65, and the velocity of rest vi sr is 86.5km.? -1 for 
the two velocity components, respectively. The data were im- 
aged with a robust weighting scheme close to uniform weight- 
ing. The synthesized beam of the line and continuum data is 
4.3" x 3.0" with a PA. of 14 degrees. The 3<x continuum rms 
is 0.24mJybeam _1 . The 3<x rms of the N2H + (l-0) data mea- 
sured from an emission-free channel with a spectral resolution of 
0.2kms _1 is 21 mJybeatrT 1 . Data calibration and imaging was 
performed with the GILD AS software package 1 . 



2.3. Nobeyama 45 m telescope 

The N2H + data has been observed using the BEARS receiver 
at the NRO 45 m telescope in Nobeyama, Japan, in April 2010. 
The receiver has been tuned to 93.17346 MHz, covering the full 
hyperfine structure of the N2H + (1 - 0) transition with average 
system temperatures of r sys ~ 200 K. At this frequency the tele- 
scope beam is 18.2" and the observing mode provides a spectral 
resolution of 0.2 km s -1 at a bandwidth of 32 MHz. The software 
package nostar (Sawada et al., 2008) was used for the data re- 
duction, sampling the data to a pixel size of 7" with a spheroidal 
convolution and smoothing the spectral resolution to 0.5 km s -1 . 
The rms of the final map is ~ 0.1 1 K per 0.5 km s -1 wide chan- 
nel. Strong winds hampered the pointing for part of the observa- 
tions, contributing to spatial uncertainties. 



http://www.iram.fr/IRAMFR/GILDAS 
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Table 1. Core parameters from the PdBI 3.2mm data 
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Notes: MRN1977 Mathis et al. (1977), OH94 Ossenkopf & Henning 1994, vl and v2 correspond to the fitted peak velocities of the two 
components in Fig. 6. 

2.4. Merging of PdBI and Nobeyama data 

To combine the N2H + (1 - 0) from the Nobeyama 45 m with the 
PdBI data, the Nobeyama data were converted into units of Jy, 
and then the combination was done within the GILDAS soft- 
ware package by the task uvshort. The synthesized beam of 
the combined dataset is 5.1" x 3.2" with a P.A. of 8 degrees. 
The 3<x rms of the merged N2H + (l-0) data measured from an 
emission-free channel with a spectral resolution of 0.2 km s~' is 
54mJybearrr 1 . 



2.5. APEX and the IRAM 30m continuum data 

The 1 .2 mm continuum data were first presented in Beuther et al. 
(2002) and the APEX 870 //m data are part of the ATLASGAL 
survey of the Galactic plane (Schuller et al., 2009). Beam sizes 
of the two datasets are 10.5" and 19.2", respectively. The l<x 
rms values are 12 and 63 mJy beam -1 , respectively. 



3. Results 

3.1. Continuum emission 

Figure 1 presents an overview of the far-infrared to submm 
continuum data obtained toward the whole cloud complex with 
Herschel and APEX. While the IRAS source and an additional 
far-infrared source in the north are bright at all presented wave- 
lengths, our target region IRDC 18310-4 shows considerable 
structure changes between 70 /mi and 870 /vm. While it is a 
strong emission source at long wavelengths, the source is only 
depicted as an absorption shadow at 70 and 100/im. While 
we cannot exclude that some ongoing low-mass star formation 
may exist below our detection threshold, non-detections at far- 
infrared wavelengths are among the best indicators for genuine 
high-mass starless gas clumps. 




18 33 40 ; 
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Table 2. Projected nearest neighbor separations 



(") (AU) 



mml-mm2 2.9 
mml-mm3 7.1 
mml-mm4 9.5 



14200 
34900 
46500 



The longest wavelengths 1.2 mm continuum data with a spa- 
tial resolution of 10.5" first presented in Beuther et al. (2002) 
are best suited to re-estimate the gas mass and column den- 
sities of IRDC 18310-4. Peak and integrated 1.2 mm fluxes 
are 132mJybeairT 1 and 717 mJy for that region, respectively. 
Assuming optically thin thermal emission from dust, a gas-to- 
dust mass ratio of 186 (Jenkins, 2004; Draine et al., 2007) and a 



Fig. 2. Zoom into the central IRDC 18310-4. The color-scale shows 
the 70yum extinction feature, and the contours present the 3.2 mm 
PdBI continuum emission starting at 3tr and continuing in 2cr steps 
(lcr = 0.08m T ybeanT'). The mm-labels, a synthesized beam and a 
scale-bar are shown as well. The phase reference center is close to mml 
at R.A. (J2000.0) 18:33:39.468 and Dec. (J2000.0) -08:21:07.65, and 
the primary beam of the observations with ~ 54" is larger than the 
shown image size. 



temperature of 18 K derived from NH3 observations (Sridharan 
et al., 2005) as well as the Herschel continuum data (see be- 
low), we can estimate masses and column densities for differ- 
ent dust properties. Using standard ISM dust properties (Mathis 
et al. 1977, MRN1977, (K L2mm ~ 0.4 cmV 1 )), the clump mass 
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and peak column density are ~16OOM and 2.5 x 10 23 cn"f 2 , 
respectively. If we use instead the dust properties discussed 
in Ossenkopf & Henning (1994) for thin ice mantles at den- 
sities of 10 5 cirT 3 (OH94, Ki.imm ~ 0.8cm 2 g _1 ), the corre- 
sponding values for mass and column density are ~800 M and 



1 .3 x 10 cm" , respectively. These two estimates can be consid- 
ered as brackets around the real gas masses and column densities 
of that region. 



Color: Tdust (K) 
Contour: 70um absor 




40s 39s 

RA (J2000) 
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Fig. 3. The color-scale presents the dust temperature map derived from 
the Herschel far-infrared continuum data. Contours outline the structure 
of the 70 fim extinction, and the stars mark the positions of the four mm- 
sources from Figure 2. 



Assuming a standard initial mass function (e.g., Kroupa 
2001) and a relatively high star formation efficiency of ~ 30% 
which should be appropriate for such dense gas clumps (e.g., 
Alves et al. 2007), one needs approximately 1000 M initial 
gas mass to form a cluster with at least one 20 M star in the 
end. Given the mass ranges estimated above, IRDC 18310-4 
should be capable of high-mass star formation. Considering col- 
umn density thresholds, Tackenberg et al. (2012) recently esti- 
mated that the initial gas column density for the Orion Nebula 
Cluster should have been approximately 1.8 x 10 23 cirT 2 , which 
is roughly agrees with the theoretical column density threshold 
from Krumholz & McKee (2008). Comparing this to the peak 
column density derived from the 1.2 mm data above, it also in- 
dicates that this region may form high-mass stars in the future. 
CommerQon et al. (201 1) recently pointed out that not only high 



column density and heating can prevent fragmentation but that 
the magnetic field may be as important. 

Figure 2 presents a zoom toward our central target region 
IRDC 18310-4. With the PdBI 3.2 mm continuum data, we can 
identify four sub-sources (mml to mm4, in the following also 
called cores, in contrast to the larger-scale environment labeled 
clump) above a 5<x threshold. While we detect a few more 3<x 
features, for the following analysis we only consider the 5<x de- 
tections mml to mm4. The cores are marginally resolved, and 
the source separation between mml to mm3 is done "by eye" at 
the lowest flux decrements between the cores. Outer edges are 
chosen at the 3<x level. Peak and integrated fluxes for these four 
cores are presented in Table 1. Using the dust emission again 
to estimate gas masses and column densities with the same as- 
sumptions outlined above (MRN1977, K3.2mm ~ 0.06 cm 2 g -1 ), 
OH94, Ki2mm ~ 0.17cm 2 g _1 ), for these four cores we find 
masses between 18 and 107 M and peak column densities be- 
tween 1.1 x 10 23 and 6.1 x 10 23 cm -2 , respectively (Table 1). 
While we obviously filter out a significant fraction of the flux 
with the interferometer (~ 73%), this region still hosts relatively 
massive cores within a high-mass gas clump reservoir. 

While the locations of mml, mm3 and mm4 are largely as- 
sociated with the lOfim extinction feature (although mml and 
mm3 are both shifted a bit to the edge), it is surprising that mm2 
lies at the edge of the dark extinction feature at a narrowing of 
the extinction. Although we cannot identify a clear and isolated 
70 pm or shorter wavelength point source at this position, the 
slight increase of 70 /mi emission toward the mm2 peak is con- 
sistent with a very young embedded object that just started to 
heat up the environment. In this picture, it may be that mm2 is in 
a little more advanced evolutionary stage compared to the other 
mm peak positions. 

Another characteristic of the fragmentation of the gas clump 
is the projected nearest neighbor separation which is listed in 
Table 2. Since the beam is elongated in approximately in north- 
south direction, we can resolve in east-west direction a separa- 
tion between mml and mm2 of about ~ 2.9", which corresponds 
at the given distance of 4.9 kpc to a projected linear separation 
of of ~ 14200 AU. The projected nearest neighbor separations to 
the other cores is even larger, between ~ 35000 and ~ 47000 AU 
(Table 2). In comparison to these values, the approximate Jeans- 
length where gas clumps should fragment into sub-structures at 
approximate temperatures and densities of 18 K and 10 5 cirT 3 is 
~ 17000 AU. Hence, we see fragmentation above and below this 
approximate Jeans-length which is reasonable since gas clumps 
have a density structure and therefore not one unique fragmen- 
tation scale. 

We also use the Herschel data to derive the temperature 
structure of the region fitting the SEDs pixel by pixel (see also 
Stutz et al. 2010 for details of the fitting). We excluded the 
500 jum data to achieve a higher spatial resolution for the temper- 
ature map of 24.9", corresponding to the beam of the Herschel 
350//m data. Figure 3 presents the results and we see a clear 
drop of the far-infrared dust temperature from the outside to the 
center of our region where temperatures around 18 K are mea- 
sured. This is very similar to the gas temperature derived toward 
that region by NH3 observations (Sridharan et al., 2005), and it 
is further evidence that no internal heating source exists so far. 



3.2. Spectral line emission 

To study the kinematics of the gas, Figure 4 presents the 
N 2 H + (1 - 0) data (critical density at 20 K is 1.6 x 10 5 cirT 3 ) of 
the region observed by Tackenberg et al. (in prep.). Although the 
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Fig. 4. The bottom panel shows in color the integrated N2H + (1 - 0) 
emission (the whole hyperfine structure between 76 and 94kms~') ob- 
served with the Nobeyama 45 m telescope toward the entire complex 
containing IRAS 18310-0825 in the south-east and IRDC 18310-4 in 
the north-west. The contours show the 1.2 mm continuum data pre- 
sented first by Beuther et al. (2002). The contouring is from 10 to 90% 
of the peak emission of 400 mjy beam 1 . The top panel shows a zoom 
into IRDC 18310-4 where the color scale presents a 1st moment map 
(intensity-weighted peak velocity) of the main hyperfine component 
complex between 84 and 90kms~'. The black contours are the same 
1 .2 mm map as in the bottom panel, and the white contours outline the 
70/jm extinction feature. The stars mark the positions of the four PdBI 
3.2mm sources. 



N2H + emission toward IRDC 18310-4 is a bit weaker than to- 
ward the HMPO IRAS 18310-0825, we nevertheless see a clear 
N2H + peak toward our target region. Zooming into IRDC 18310- 
4, the top panel of Figure 4 shows a 1st moment map (intensity- 
weighted peak velocity) extracted from these data using only the 
main hyperfine structure peak (see also Fig. 5). This figure shows 
on the one hand that the overall velocity structure of the region 
is relatively uniform with no pronounced velocity structure and 



only weak gradients across the whole region (< lkms" 1 over 
~1 pc extent). 

Figure 5 shows a spectrum extracted toward the peak posi- 
tion. Observing the whole N2H + (1 - 0) hyperfine structure si- 
multaneously allows us to fit the optical depth and by that derive 
accurate spectral line parameters. For the single-dish spectrum 
toward the peak position, we find a line-width of 2.7 km s -1 at a 
velocity of 86. 1 km s . The fitted peak velocity is very close to 
that measured previously in NH3 (86.5 km s , Sridharan et al. 
2005). Following MacLaren et al. (1988), we can estimate a 
virial mass assuming the FWHM/2 of the Nobeyama beam as 
the radius of the core. Employing possible density structures of 
1/r and 1/r 2 as border conditions, the virial mass of this region 
is between 305 and 202 M Q , respectively. Compared to the to- 
tal gas masses derived from the single-dish dust continuum data 
described in the previous section (section 3.1), the virial mass is 
considerably lower (between factors 2.6 and 7.7, depending on 
the assumptions for the dust properties and density structure). 
This implies that the whole region is gravitationally unstable on 
these scales and likely already undergoing collapse. 



Fig. 5. Compilation of N2H + (1 - 0) Nobeyama 45 m single-dish spec- 
trum toward IRDC 18310 (black histogram with green fit) and PdBI 
spectrum toward mml (red histogram). The PdBI data are multiplied 
by 5 for clarity. Fit results to the whole N 2 H + (1 - 0) hyperfine structure 
are given in the panel. 



Since our PdBI observations are not only sensitive to the 
continuum but also to the N2H + (1 - 0) emission, we have kine- 
matic information at high spatial resolution as well. Figure 5 
also shows an N2H + (1 - 0) spectrum extracted toward the peak 
of mml just from the PdBI data (scaled up by a factor 5), and 
one clearly sees kinematic differences between the large and the 
small scale. To account for the differences and to complement 
the PdBI data with the missing flux information, we merged the 
PdBI and Nobeyama data into one single data-cube (see Section 
2.4). 

Figure 6 presents the spectra extracted from the merged data- 
cube toward the mm peak positions. All these spectra show that 
we have at least two separated velocity components in the re- 
gion, one approximately peaking at ~86.5kms~ 1 and the other 
at ~89.3kms~ I (these values are averages from the full hyper- 
fine structure fits to mml to mm3, see Fig. 6). 



Fig. 6. Spectra and N 2 H + (1 - 0) hyperfine spectral line fits extracted 
from the PdBI+Nobeyama N 2 H + (1 - 0) data toward the mm-peak posi- 
tions marked in Figure 2. Fit results and source labels are presented in 
each panel. 



The comparison of the PdBI only spectrum and the 
Nobeyama only spectrum in Figure 5 shows that the lower- 
velocity PdBI component corresponds relatively well to the 
single-dish spectrum whereas the higher-velocity PdBI compo- 
nent does not have a proper single-dish spectrum counterpart. 
This strongly indicates that the ~86.5kms _1 component is the 
main component also present in the ambient gas, whereas the 
higher velocity component at ~89.3 km s _1 appears to be mainly 
associated with a dense sub-part moving with respect to the am- 
bient gas. 

To visualize spatial differences between these two spectral 
components, we imaged the N2H + data for both velocity struc- 
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Fig. 7. Pdbl only (left panel) and merged PdBI+Nobeyama 45 m N 2 H + (1 -0) data (right panel). The color-scale shows the 70//m extinction feature, 
and the contours present the integrated N2H + (1 - 0) emission of three different velocity components as indicated in the figure. To avoid hyperfine- 
splitting overlap, these integrated maps were created from the isolated hyperfine component -8.0kms~' apart from the central component, and 
then the velocity shifted to the center to be comparable with other line measurements. The contouring is done from 15 to 95% of the peak emission 
of each dataset, respectively. The mm-labels, a synthesized beam and a scale-bar are shown as well. 



tures separately. Figure 7 shows the integrated emission derived 
for the velocity regimes [85.5,88.5] and [88.5,89.7] km s -1 , re- 
spectively. The left panel shows the PdBI- only data whereas 
the right panel presents the combined PdBI+Nobeyama 45 m 
data. To avoid any artifacts of overlapping hyperfine structure 
lines, we produced these images from only the isolated hyperfine 
structure component offset by -8.0kms _1 (F — - 1). While 
the lower-velocity component agrees well with the mid-infrared 
absorption structure encompassing mml, mm3 and mm4, the 
higher-velocity component is peaking mainly in the direction 
of mm2. It is at first sight surprising to note that the higher- 
velocity component in the PdBI-only data shows an extension to- 
ward mm4 whereas this is lost in the merged dataset. Most likely 
this is because that extension is on small scales and compara- 
bly weak that it gets lost by the merging process. Additionally, 
it is possible that the single-dish data add more noise in the 
combined dataset which make this weaker extension harder to 
detect. In addition to this, the merged dataset shows one ad- 
ditional velocity component at even higher velocities between 
89.7 and 90.6kms~ 1 that is more diffusely distributed north of 
the main IRDC 18310-4 structures. It is unclear whether this 
is just another unrelated gas component or whether that gas is 
structurely related to the red-shifted component between 88.5 
and 89.7 km s~'. Nevertheless, both main velocity components 
strongly overlap and hence both contribute to the N2H + emis- 
sion from mml and mm2 (see spectra in Fig. 6). Furthermore, 
the PdBI-only data indicate that mm4 is located at an interface 
region between both components. In contrast to this, mm3 ap- 
pears a bit different because it is located at the edge of the inte- 
grated emission of both components. 

Table 1 lists the line widths of the two components as well 
as the estimated virial masses following MacLaren et al. (1988) 
for different density profiles (p oc - and p oc 4)- The assumed 
radius of the sources is 1/2 of the FWHM which corresponds 
at the given distance to ~ 9000 AU. It is obvious that the nar- 
row higher-velocity component barely contributes to the stabil- 
ity of the cores, and that the broader low-velocity component is 
far more important to stabilize the structures. Considering the 
spread in virial masses for different density profiles as well as 
for the gas masses based on the dust continuum emission and 
the assumed different dust properties (Table 1), the cores appear 
closer to virial equilibrium. This is different to the overall virial 
to gas mass ratio based on the single-dish data for the whole re- 
gion that is significantly smaller than 1 (see above). However, 
one has to keep in mind that the virial masses are calculated 
from the merged PdBI+Nobeyama N2H + data whereas the core 
masses are calculated from the PdBI data only. If one compares 
the main component of the N2H + spectrum toward mml without 
and with the Nobeyama short spacing information (Figs. 5 & 6), 
one finds that the combined data reveal a FWHM about a fac- 
tor 1 .7 higher than the interferometer data only. Since the line- 
widths is squared in the estimate of the virial mass, the derived 
virial masses of the combined data are usually a factor of a few 
larger than with interferometer data only. This implies that the 
comparison of masses from the combined PdBI+Nobeyama line 
data with the PdBI-only continuum data is not reliable. With the 
missing flux problem, the core masses from the PdBI-only data 



are underestimated, implying that the ratio of virial to gas mass 
on the core scales should be smaller than 1 as well. Hence, the 
cores are also prone to collapse. 

Figure 8 presents the integrated N2H + emission as well as the 
first and second moment maps (intensity-weighted peak veloc- 
ities and line widths). In particular the 1st moment map shows 
that we have no strong velocity gradients across the filament but 
that the velocity peaks are dominated by the two components 
discussed above. The 2nd moment map appears as if it had a line 
width increase toward mm2, however that is only because the 
2nd moment algorithm assumes a single component, and there- 
fore the two real components appear as one broader feature in 
that direction. 



4. Discussion 

4.1. Linear scales and densities from low- to high-mass 
starless gas clumps 

How do linear scales and average densities compare between 
star-forming regions from low- to high-mass? To compare the 
different mass regimes with the high-mass starless region pre- 
sented in this paper, we selected the low-mass starless core 
B68 (e.g., Alves et al. 2001; Nielbock et al. 2012) which re- 
sembles typical other low-mass globules (e.g., Launhardt et al. 
2010), and the intermediate-mass starless star-forming region 
IRDC 19175 (Beuther & Henning, 2009). Figure 9 presents data 
for all three regions. It is interesting to note that all three re- 
gions have approximately the same linear spatial extent on the 
order of 50000 AU. However, while the linear size between the 
three regions does not vary considerably, the total gas reser- 
voirs for the three complexes are very different, ranging from 
-3.1 M Q for the low-mass core B68 (Nielbock et al., 2012), 
~87M Q for the intermediate-mass region IRDC 19175 (Beuther 
& Henning, 2009) to values between -800 and ~1600M o (de- 
pending on the dust model, see section 3.1) for the high-mass 
region IRDC 18310-4 (see Table 3 for a list of the source param- 
eters). Assuming that the three regions resemble typical starless 
gas clumps in the different mass regimes, these data imply that 
not the linear sizes are important but the mass reservoir squeezed 
within such regions, i.e., the densities. 

While the peak density of the B68 core is ~ 3.4 x 10 5 ctrT 3 
(Nielbock et al., 2012), the average density of the region is con- 
siderably lower. We estimate an average density for B68 from the 
total mass of 3.1 M Q (recent analysis of Herschel data, Nielbock 
et al. 2012) divided by a sphere with approximate radius of 
~ 150" at a given distance of ~ 150 pc. Using these numbers the 
average gas density of the low-mass starless core B68 is approx- 
imately 1.2 X 10 4 cirT 3 over a region corresponding to a sphere 
with -45000 AU diameter. 

For the intermediate-mass star-forming region IRDC 19175, 
average core densities are estimated by Beuther & Henning 
(2009) to ~ 6 x 10 5 cm~ 3 over regions corresponding to spheres 
with ~3500AU diameter. Estimating an average density from 
the total mass of 87 M over a volume corresponding to an 
effective radius of the 1.2 mm continuum map in Beuther & 
Henning (2009) results only in slightly lower values around 
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Fig. 8. N 2 H + (1 - 0) moment maps of the PdBI-only data toward IRDC 18310-4. The left, middle and right panel show in color the 0th moment 
(integrated intensity), 1st moment (intensity-weighted peak velocity) and 2nd moment (intensity-weighted linewidth), respectively. The moments 
were again extracted for the isolated hyperfine component -8.0km s~' from the line center. The 1st moment was calculated from this isolated 
hyperfine component and then shifted toward by +8.0 km s -1 to the v/ sr of the source again. The white contours outline the 70 //m extinction 
feature and the stars mark the 3.2 mm peak positions from Figure 2. The synthesized beam and a scale-bar are shown again. 

Fig. 9. Comparison between starless gas-clumps from low- via intermediate- to high-mass star formation. The 875 /jm data for the prototypical 
low-mass core B68 are taken from Nielbock et al. (2012), and the N 2 H + (l-0) data for the intermediate-mass star-forming region IRDC 19175 are 
from Beuther & Henning (2009). Each panel presents a linear scale-bar of 50000 AU and the synthesized beam, average densities are mentioned 
as well. 



~ 2.6 x 10 5 cm 3 over a region corresponding to a sphere with 
-48000 AU diameter. 

To compare these values with our high-mass star-forming re- 
gion IRDC 18310-4, we estimate an average density for the core 
mml via averaging the peak column densities listed in Table 1 
for the two different dust models, and assuming that the 3rd di- 
mension has the same size as our synthesized beam. This results 
in average densities for mml of ~ 1.5 x 10 6 ctrT 3 over a re- 
gion corresponding to a sphere with ~ 18000 AU diameter. To 
get the density over a slightly larger regions, we use the average 
column density derived in section 3.1 for the two different dust 
models from the single-dish 1.2 mm continuum data and again 
assume the 3rd dimension to have the same spatial extent as the 
single-dish beam of 10.5". This results in average densities of 
~ 2.5 x 10 5 ctrT 3 over a region corresponding to a sphere with 
~5 1 000 AU diameter. 



Table 3. Star-forming region parameters 





B68 


IRDC19175 


IRDC 183 10-4 


mass [M Q ] 


3.1 


87 


800-1600 


av. dens, [cm -3 ] 


1.2xl0 4 


2.6xl0 5 


2.5xl0 : ' 


size [AU] 


45000 


48000 


51000 


Jeans length 2 [AU] 


44000 


10000 


10000 


Jeans mass 2 [M Q ] 


1.7 


0.36 


0.37 


core dens, [cm -3 ] 


_i 


6.0x1 5 


1.5xl0 6 


Jeans length 2 [AU] 


_i 


6000 


4000 


Jeans mass 2 [M ] 


_i 


0.24 


0.15 


1 B 68 as a whole is already the core. 




2 At 15 K 









These large density differences obviously also have a strong 
impact on the fragmentation properties of the starless gas 
clumps. While the low-density low-mass core B68 does not frag- 
ment at all, the higher-density intermediate- to high-mass star- 
less clumps do fragment significantly. To better understand these 
properties, we estimate the Jeans length for the average B68 den- 
sity of 1.2xl0 4 cm -3 , the approximate average large-scale densi- 
ties of IRDC 19175 and IRDC 18310-4 of ~ 2.5 x 10 5 cm 3 , and 
the higher core densities of the two regions of ~ 6 x 10 5 cm -3 
and ~ 1.5 x 10 6 ctrT 3 , respectively. Assuming average temper- 
atures of 15 K for all these starless regions, the corresponding 
Jeans-length are -44000 AU (B68), ~ 10000 AU (average large- 
scale densities of IRDC 19175 & IRDC 18310-4), -6000 AU 
(average core density of IRDC 19175) and ~4000 AU (average 
core density of IRDC 18310-4), respectively. This relatively sim- 
ple isothermal Jeans-analysis already explains the observations 
well: The large Jeans-length for B68 is reflected in the fact 
that this low-mass starless core does not fragment on the given 
scales (Fig. 9). Furthermore, as already discussed in Beuther 
& Henning (2009), the different Jeans-lengths for B68 (located 



close to the Pipe nebula) and IRDC 19175 correspond well to 
the observed core separations in the low-mass Pipe nebula and 
the intermediate-mass IRDC 19175 region, respectively. For the 
high-mass region IRDC 18310-4, we are still spatial resolution 
limited and only observe the fragmentation on the scales of our 
observed resolution elements (see also Table 2). In future higher- 
spatial-resolution observations, we expect to resolve this region 
into even more and closer spaced cores. 

Calculating the Jeans-mass instead of the Jeans-length for 
the average densities of B68, IRDC 19175 and IRDC 18310-4, 
we find 1.7, 0.36 and 0.37 M Q , respectively. The Jeans-masses 
for the higher core densities are correspondingly smaller (Table 
3). For the low-mass region B68 that value is on the order of 
the whole mass (3.1 M Q ), and also for the intermediate-mass re- 
gion IRDC 19175, the estimate masses of the individual cores 
are mostly below 1M . Hence Jeans-masses and core masses 
approximately agree in these two cases. In contrast to that, the 
Jeans-mass in IRDC 183 10-4 is significantly lower than the mass 
estimates we derived from the dust continuum data (Table 1). 
Therefore, similar as discussed in the previous paragraph, for 
the high-mass region IRDC 18310-4, we are still spatial resolu- 
tion limited, and it is highly likely that the here presented cores 
mml to mm4 will fragment even further when observed at higher 
spatial resolution. 

Assuming that these example regions resemble starless gas 
clumps over the whole mass regime, one may conclude that 
gravitational instability in the classical Jeans-analysis explains 
the fragmentation properties of starless gas clumps from low- to 
high-mass reasonably well. The major difference between the re- 
gions appears to be their densities and mass reservoirs squeezed 
into the respective spatial scales. Similar Jeans fragmentation 
scales have also been observed toward more evolved regions, 
e.g., the young cluster NGC2264 (Teixeira et al., 2006). In con- 
trast to the larger-scale cloud fragmentation where turbulence 
plays a major role (e.g., Padoan & Nordlund 2002; Mac Low 
& Klessen 2004), these data indicate that the fragmentation on 
smaller clump and core scales appears to be dominated by grav- 
ity. 

4.2. Kinematic signatures of a dynamic cloud collapse 

Independent of the cloud formation models - do clouds evolve 
slower via ambipolar diffusion processes (e.g., Mouschovias 
et al. 2006; McKee & Ostriker 2007) or faster via more dy- 
namic gas flows (e.g., Vazquez-Semadeni et al. 2006; Heitsch 
et al. 2008; Banerjee et al. 2009) - the gas has to show dynamic 
signatures of contraction and/or infall over large spatial scales. 
During molecular cloud formation the diffuse HI gas converts 
to the denser molecular phase. Within molecular clouds, fila- 
mentary sub-structures are likely important, and recent data in- 
dicate that hierarchical, often also filamentary structures feed the 
central and densest filaments or gas clumps (e.g., Myers 2009; 
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Schneider et al. 2010; Schmalzl et al. 2010; Hennemann et al. 
2012). These central cluster-forming clumps have then to col- 
lapse and form the multiple and clustered stellar associations 
in the end. Statistical signatures of infall from single-dish aver- 
age line profiles exist. However, they are often ambiguous, and 
one finds signatures for infall as well outflowing motions (e.g., 
Fuller et al. 2005). More coherent infall motions on larger scales 
have only rarely been reported (e.g., Motte et al. 2005; Schneider 
et al. 2010). Going from the cluster- forming scales down to cores 
forming individual or multiple, bound objects, dynamic studies 
in the past have mostly concentrated on more evolved evolu- 
tionary stages like high-mass protostellar objects (HMPOs) or 
ultracompact Hn regions (e.g., Keto 2002; Sollins et al. 2005; 
Sandell et al. 2009; Qiu et al. 2011; Goddi et al. 2011; Beuther 
et al. 2012a). 

While the above infall studies dealt with the average infall 
motions of the entire star-forming gas clumps, only recently 
spatially resolved dynamical studies of sub-structures and cores 
within the natal gas became possible (e.g., Beuther & Henning 
2009; Zhang & Wang 2011; Wang et al. 2011; Csengeri et al. 
2011a,b). In contrast to the above discussed single-dish studies 
where often infall signatures are derived from line profiles of 
optically thick lines, Ragan et al. (2012b) also find signatures 
of collapse and fragmentation via high-spatial-resolution NH3 
observations of IRDCs. Furthermore, the studies by Csengeri 
et al. (2011a,b) revealed multiple velocity components in opti- 
cally thin gas tracers like H 13 CO + or N2H + that were interpreted 
as signatures of shear motions possibly caused by converging 
gas flows. 

The spectra we found (Figs. 5 and 6) also show multiple ve- 
locity components in the dense gas tracer N2H + . Fig. 7 even 
identifies that the different velocity components are dominated 
by different mm cores within the large-scale gas clump. In addi- 
tion to the work by Csengeri et al. (2011a,b), multiple velocity 
components in optically thin dense gas tracers were also found 
by Beuther & Henning (2009) as well as in recent work by Ragan 
et al. (in prep.) with the PdBI toward the Gl 1. 1 1 filament, and by 
Bihr et al. (subm.) with VLA NH3 observations toward IRDCs. 
Combining our new data with the results recently found by the 
other groups, based on high-spatial resolution interferometer ob- 
servations there emerges a regular signature of multiple veloc- 
ity components in optically thin dense gas tracers. Typical sep- 
arations line widths of these multiple velocity components are 
around ~2 km s -1 and ~1 km s -1 , respectively. 



Fig. 10. The N2H + (1-0) isolated hyperfine line profiles viewed through 
various lines of sight through a massive star forming region. The central 
image with a side-length of 0.8 pc shows the column density of the re- 
gion and the surrounding line profiles are calculated for a 0.09 pc beam 
centered directly on the embedded core and viewed along the direction 
each panel touches the central image. 



Confronting these observational results with theory, Rowan 
Smith has recently started to simulate the spectral line emission 
from dynamically collapsing clumps of a large-scale collapsing 
gas cloud modeled initially by Smith et al. (2009) and Bonnell 
et al. (2011). First results about cores associated with low-mass 
star formation were presented in Smith et al. (2012), and the sim- 
ulations for high-mass collapsing gas clumps are currently con- 
ducted (Smith et al. subm.). Figure 10 shows an example of non- 
Gaussian optically thin line profiles produced from a simulation 



of a massive star-forming gas clump with 370 M Q (Smith et al. 
2009, Smith et al. subm.). The N 2 H + (1-0) line profiles were 
calculated using an assumed abundance ratio of 10~ 10 relative to 
molecular hydrogen (Aikawa et al., 2005). The radiative trans- 
fer code RADMC-3D (Dullemond, 2012) was used to calculate 
the line emission given the temperatures and densities from the 
original simulation. (For a full description of our method see 
Smith et al. 2012). Here, only the isolated hyperfine component 
F - 0- 1 (see Fig. 6) offset by ~8.0km s from the central com- 
ponent is shown. The line profiles are synthesized over a beam 
of FWHM 0.09 pc and calculated for a variety of viewing angles. 
No artificial noise is added to the data. 

Clearly, these simulations find similar signatures of multiple 
N 2 H + components within the dynamically collapsing gas clump. 
Typical peak separations and line widths are similar to what is 
found in the observational data. The multiple components seen 
in Figure 10 arise in the simulations due to the clumpy nature 
of the high-mass star-forming region. Across the region there is 
a large scale supersonic collapse motion running from around 
+2kms _1 on one side to -2kms _1 on the other. Dense clumps 
of gas at different positions along this velocity gradient all con- 
tribute to the total emission and produce peaks in the line profile 
at their individual velocities. A full analysis of line profiles from 
massive star formations will be presented in Smith et. al.(subm.). 

This interpretation is also consistent with our finding that the 
ratio of viral mass to total gas mass is significantly smaller than 
1 for the entire region, implying that the large-scale gas clump is 
prone to collapse, similar to the findings of Ragan et al. (2012b) 
for their sample of IRDCs. Since the virial mass to core mass 
ratio is smaller than 1 on the core scales as well, the region is 
consistent with a globally collapsing gas clumps on almost all 
scales. 

We can also estimate the dynamical crossing time as well 
as the free-fall time of the region. Using the typical velocity 
difference between the two gas components of ~ 2.8kms~' as 
the velocity dispersion of the clumps, and the approximate ex- 
tent of the 70//m absorption feature of ~ 30" (corresponding to 
147000 AU at the given distance), the dynamical crossing time 
for IRDC 18310-4 is ~ 2.5 x 10 5 yrs. In comparison to that, es- 
timating the free-fall time for a 1000M o gas clump of with a 
radius of half the previously used extend, we get an approxi- 
mate free-fall time-scale of ~ 1.2 x 10 5 yrs. Although there is a 
difference of about a factor 2 between the two estimates, these 
are just approximate values and they agree order-of-magnitude 
wise. This again implies that a dynamical collapse of the whole 
gas clump is a feasible model for this region. 

Combining our observational results with the modeling of 
line profiles for dynamically collapsing gas clumps, we identify 
signatures that even this very young, still starless high-mass gas 
clump is actively collapsing and thus at the verge of star forma- 
tion. The fact that similar signatures were recently found toward 
other regions as well, indicates that multiple velocity compo- 
nents of optically thin gas tracers are an excellent tool to iden- 
tify and study the earliest stages of dynamical infall in high-mass 
star formation. 

5. Conclusions 

Since the evolutionary time-scale for high-mass starless gas 
clumps is extremely short (e.g., Motte et al. 2007; Russeil et al. 
2010; Tackenberg et al. 2012), so far they had barely been stud- 
ied in great depth at high spatial resolution (some notable recent 
papers in that field are Csengeri et al. 201 la; Ragan et al. 2012b). 
To better understand the fragmentation and dynamical properties 
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of the most massive gas clumps during their earliest evolutionary 
stages, here we present an interferometric PdBI kinematic and 
continuum study of a prototypical high-mass starless gas clump 
that remains far-infrared dark up to 100//m wavelengths. 

The gas reservoir between ~800 and ~1600M o (depending 
on the assumed dust properties) splits up at a spatial resolution 
of 4.3" x 3.0" (linear scales 21000 x 15000 AU) into 4 cores. A 
comparison of this high-mass starless region with typical low- 
mass (B68) and intermediate-mass (IRDC 19175) star-forming 
forming regions reveals that the overall spatial size of the low- 
mass to high-mass star-forming regions does not vary signifi- 
cantly. While the linear extent of the regions is similar, the mass 
reservoir squeezed into this spatial size on the order of 0.25 pc, 
and hence the densities do vary by about 2 orders of magnitude. 
A classical gravitational instability Jeans analysis is capable to 
explain the fragmentation properties of the gas clumps over the 
whole mass range presented. 

In addition to the fragmentation properties, we find multi- 
ple spectral velocity components toward the resolved gas cores. 
Recent radiative transfer hydrodynamic simulations of dynam- 
ical collapsing high-mass gas clumps are consistent with these 
multiple velocity components. This also agrees with a ratio be- 
tween estimated viral and total mass of the whole region <1, 

In summary, although this massive gas clump does not have 
any embedded protostellar source down to our Herschel far- 
infrared detection limits, the fragmentation and dynamical prop- 
erties of the gas and dust are consistent with early collapse mo- 
tion and clustered star formation. 
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